Ferromagnetic Resonance in Spinor Dipolar Bose-Einstein Condensates 
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We used the Gross-Pitaevskii equations to investigate ferromagnetic resonance in spin-1 Bose- 
Einstein condensates with a magnetic dipole-dipole interaction. By introducing the dipole interac- 
tion, we obtained equations similar to the Kittel equations used to represent ferromagnetic resonance 
in condensed matter physics. These equations indicated that the ferromagnetic resonance originated 
from dipolar interaction, and that the resonance frequency depended upon the shape of the conden- 
sate. Furthermore, spin currents driven by spin diffusions are characteristic of this system. 

PACS numbers: 03.75.Mn, 03.75.Nt 



I. INTRODUCTION 

Magnetic resonance (MR) as a physical concept has 
been applied in various fields, enabling physical, chem- 
ical, and medical experiments to obtain information on 
nuclear spin and electron spin systems. The concept has 
also provided valuable information to help understand 
the unknown structures of many condensed matter sys- 
tems [l|. 

The use of MR in the study of ferromagnets, e.g. 
Nickel, Cobalt, and Iron, began in the 1940s. Grif- 
fiths observed that the Lande's <?-factor of electrons in 
ferromagnets was far from the well known value, 2 
In order to understand these anomalous results, Kittel 
theoretically introduced a demagnetizing field into the 
equation representing the motion of the magnetization 
M = (M x , My, M z ), obtaining an equation valid in an 
external magnetic field HqZ, with M z q = Hq/N z and de- 
magnetizing fields fsj , thereby obtaining the Kittel equa- 
tion, 



dM 



7n[M x H] 



(1) 



Here, 7„ is the nuclear gyromagnetic ratio, and H = 
(—N X M X , —N y My, Hq — N Z M Z ) is given by the demagne- 
tizing factors iVj. By linearizing the magnetization M = 
Mo+^M from the stationary magnetization Mo = M z qz, 
Kittel obtained a precession of the magnetization and a 
precessing frequency, i.e. resonance frequency, 

^ 2 = jl{H +(N y -N z )M z0 }{H +(N x -N z )M z0 }, (2) 

which explained the anomalous <?-factor. Furthermore, 
he found that the resonance frequency depends on the 
shape of a ferromagnet because N t depends on the shape 
Thus, ferromagnetic resonance (FMR) was estab- 
lished, and the work enabled numerous additional studies 
1- 

MR also plays an important role in quantum conden- 
sate systems. In superfluid 3 He, the dynamics of the spin 
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vector and the d-vector are represented by the Leggett 
equation, which couples these vectors through magnetic 
dipole-dipole interactions Q. The equation also shows 
not only an MR typical of condensed matter, but also a 
new MR that cannot be described using the equations of 
motion for general paramagnets and ferromagnets. This 
MR was used to find A and B phases [fj[ ■ Parallel ring- 
ing, which is an oscillation of longitudinal magnetization, 
was also observed Q- 

Since the discovery of atomic Bose-Einstein conden- 
sates (BECs) Hi], BECs have been studied in optics and 
atomic and condensed matter physics. We have intro- 
duced MR into BECs to realize magnetic resonance imag- 
ing, a popular method of nondestructive testing. Spinor 
BECs are expected to be suitable for MR, since they 
have not only internal degrees of freedom but also mag- 
netic properties. In particular, we are interested in mag- 
netic dipole-dipole interactions (MDDI) in spinor BECs, 
which have been actively studied. The interaction be- 
tween spins has a characteristic symmetry of rotation 
and spin, which is expected to result in a new quan- 
tum phase and Einstein-de Haas effects [l3| . Ex- 
perimentally, Griesmaier et al. realized spinor dipolar 
condensates using 52 Cr atoms, which have a larger mag- 
netic moment than alkali atoms [14]. The shape of the 
condensates clearly represented the anisotropy of the in- 
teraction [THHH. Thus, MDDI has opened new areas of 
spinor condensate research. 

As an introduction to MR in BECs, we numerically 
studied spin echo in dipolar BECs with spin-1 [19]. The 
spin echo is a typical phenomenon of MR, discovered by 
Hahn [13] and developed by Carr and Purcell [H[ . Previ- 
ously, we calculated the transition from Rabi oscillations 
to internal Josephson oscillations in spinor condensates 
[2fj| . In this paper, we consider MDDI in spin-1 BECs, 
examining FMR by analyzing the Gross-Pitaevskii (GP) 
equations. 

In section [H] we derive Kittel-like equations from the 
GP equations, and analyze them. In section IIII1 using 
a single-mode approximation, we derive Kittel equations 
from the Kittel-like equations. The MDDI of the Kittel 
equations is considered as the origin of the demagnetizing 
field, which is phenomenologically introduced in Eq. ([1} . 
In section IIV1 we numerically solve the GP equations, 
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obtaining resonance frequencies that depend upon the 
shape of the condensates, and spin currents driven by 
spin diffusion which is given by the MDDI. Finally, Sec. 
IVlis devoted to our conclusions. 



II. FORMULATION 

In this section, we derive the equations of motion for 
spins from the spin-1 GP equations with an external mag- 
netic field and an MDDI 
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Here, V is the trapping potential, \x is the chemical 
potential, and the total density n — J2i n i i s gi yen 
by rii = |0i| 2 . The external magnetic field is H = 
(H x , H y , H z ), and the components F l a o of the spin ma- 
trices Fi are for spin-1. The interaction parameters are 
co = (.9o + 25 2 )/3 and c 2 = (ff2-go)/3 for gi = Airh 2 a l /M 
represented by s-wave scattering lengths a, . The dipolar 
coefficient is Cdd = MoSe^s/"^, and the unit vector is 
e = (e x , e y , e z ) — (x — x', y — y', z — z')/\v — r'\. 

Under the homogeneous magnetic field H = Hi, the 
equations can be rewritten as, 
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dt V 2M 
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(4a) 



(4b) 
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+ c 2 {(n_i + n - ni)V>-i + V'iV'o} + D -i 
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These dipolar terms are represented as, 
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with the integrations c?± = d x ± idj, and d z given by, 

^ = ^j^y^{ l ~^H^- ( 5 ) 



The spin density vectors are defined as, 

= A{V,*(^ 1 +0_ 1 )+^(V'i+C 1 )}, (6a) 
F, = &F y * 

= -/={0o(V>x -V>-x) - V'oCV'i -V'-i)}) ( 6b ) 

F z = *t^ = R (|^ 1 |2_|^_ 1 |2j - (6c) 

Here, ^ = (-01, 0Oj 0-i) T is the spinor wave function. 

Differentiating Eq. ([6]) with respect to time and utiliz- 
ing Eq. (j4j, we can obtain the Kittel-like equation, 



<9F 

— = K + 7e [FxH cff ] 
at 



(7) 



with the gyromagnetic ratio 7e = g^is/h of an electron. 
The first term K = (K x , K y , K z ) becomes, 

h 1 



K. 
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The effective magnetic fields H c ff = H + Hdd = 
(H* s , H^ s , H* s ) consist of the external magnetic field and 
the dipolar field Hdd, given by, 

ttx _ C dd 7 

-"eff — 
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Note that Eq. (O does not depend on spin exchange 
interaction, which refers to the second term with c 2 in 
Eq. ([3]). Generally, the interaction affects a spin through 
the effective magnetic fields of the other spins. However, 
exchange interaction does not appear in H e g. Therefore, 
the isotropic exchange interaction does not affect MR in 
these condensates. 

We can redefine Eq. ([7]) as, 

dF h 

ST = 7777TV 2 F fc -V-j fc +7e[FxH off ] fc , (8) 



dt 2Mi 



where, 



J.c 



;(VSV(^ + 0-l) + (04 + 0* i)V0 O ), 
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h - 7f^^oV(0i -0- 1 )-(0 1 *-01 1 )V0 o ), 
3* - ^(VWx-0*iV0_i). 
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The equation of motion @ for spins describes the prop- 
erties of spin dynamics in a ferromagnetic fluid. The 
first, second, and third terms of Eq. ([8]) represent spin 
diffusion, spin current, and spin precession around H e ff, 
respectively. 

Comparing Eq. © with Eq. (JTJ> , we noticed several 
differences. First, Eq. ([5]) was directly derived from 
the GP equations, whereas Eq. (TTJ) is a phenomenologi- 
cal equation of magnetization. The spin density vectors 
in Eq. ([5]) are microscopically affected by other spins 
through the dipolar fields in the effective magnetic fields. 
On the other hand, the magnetization in Eq. (fTJ) is af- 
fected by demagnetizing fields originating from macro- 
scopically polarized magnetization in the condensed mat- 
ter. Namely, Eq. ([8} can describe the macroscopic de- 
magnetizing field resulting from the microscopic dipolar 
field. This is a very important difference between these 
equations. 

We initially investigated the physics of the first and 
second terms of Eq. © . To simplify the discussion, we 
considered the equation under the condition H e gp = 0. 
Thus, we derived the continuity equations, 

BF 

^+V-J i = 0, (9) 
where J& = jfc — H/(2Mi)'VFk is an effective current term, 



:{^o*V(Vi + iM) + (Vi* + Ci)V^o 



2 V / 2M 

,r « - V-i) - (Vi " ^li)VVo 



2V2M 

ih 2 



(10a) 



(10b) 



(10c) 

Equation ([9]) can also be rewritten as, 



- J FidV = J V • J<dV = J Ji • ndS, 

by using the volume integral and the surface integral, 
whose unit vector n is vertical to the surface for Stokes' 
theorem. The equation indicates that the expectation 
value of the spin matrix (Fi) — J dVFi in the volume V 
is conserved for the spin probability flux Jj leaving and 
entering the surface. 

Under H c g 7^ 0, the Kittel-like equation can be re- 
duced to the following equation, 

BF 

■^+V.J, = [FxH eff ] 1 , (ff) 

where the right side of the equation breaks the conserva- 
tion law of spin density. Therefore, the Kittel-like equa- 
tions have two dynamics: spin precessions with frequency 
given by the effective magnetic field and spin currents 
without spin conservation. The spin currents of the sys- 
tem will be discussed in Sec. IIVBI 



III. FMR UNDER SINGLE-MODE 
APPROXIMATION 

In order to study the basic properties of the second 
term in Eq. ([7|). we introduced the single- mode approxi- 
mation, 



V>i(r,i) = VN^{t)<f>(r)exp 



i/it 



(12) 



where <j> satisfies the eigenvalue equation (-fi 2 V 2 /2M + 
V + c n)<p = pup with the relation J dr\<fi\ 2 = 1. The 
approximation is effective when the shapes of the con- 
densates are determined by the spin-independent terms, 
namely |cn) > |oa] [2l|. For 87 Rb and 23 Na, the re- 
lation is satisfied. Under this approximation, the first 
term of Eq. (|7|) vanishes, and we obtain the Kittel equa- 
tion for the spatially independent spin density vector 

S = (S X , Sy, S Z ), 



where, 



f =7e[SxH-] 



s v = 4{£&-£-i)-£o(£-ei)}, 
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and the effective magnetic field 
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(14) 

Equation ([13")) also indicates that the spin vector S pro- 
cesses around H^g IA . The precession frequency reveals 
the characteristic dynamics. Next, we consider a small 
deviation SS = (8S X , SS y ,SS z ) around the stationary so- 
lution, S = S z with So = H /Nl d , of Eq. ([15]). namely 
S = So + <5S. Introducing this representation into Eq. 
(|13p and linearizing the equation, we derived the follow- 
ing equations, 



d 
dt 
d 

d_ 

lit 



SS X = le {H + (N y dd - m d )S }SS y , 
6S V = - le {H+(NS d -NZ d )S }6S x 
tiS-_ = 0, 



which give the resonance frequency, 
u? = 7 2 e {H+(Nl d -Ni d )S }{H+(Ny d -Ni d )S } (15) 

The spin precesses with the resonance frequency lj, which 
depends on the dipolar terms N dd . 
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Here, we consider the single particle density distri- 
bution |^(r)| 2 oc e~( x +v +XzZ )/ a , where A z is the as- 
pect ratio, and discuss simple situations. For the spher- 
ical case of \ z = 1, the integration flH} results in 

^dd = ^dd = JVdd. g ivin § w = 1* H - The dipolar 
fields are canceled because of the isotropy, so that the 
spin precesses with Larmor frequency. For the circular 
plane (infinite cylinder) case of A 2 = oo (0), we obtain 
u = le {H - (iV| d - iV| d )5 } for iV| d = AT« d . 

In this representation, it seems that the microscopic 
dipolar fields, Eq. (fT4|) . act as a macroscopic demagne- 
tizing field to compare Eq. © with (|T5l) . We believe that 
the origin of the demagnetizing field is an MDDI. If the 
above discussion is correct, the dipolar coefficients iV dd 
should depend on the shape of the condensates. However, 
the single-mode approximation in spinor dipolar BECs 
is not effective in large-aspect-ratio condensates, as dis- 
cussed by Yi and Pu |22j. Therefore, we must consider 
the spin dynamics beyond the approximation. 

IV. FMR FOR NUMERICAL CALCULATION 
A. Precession dependence on the aspect ratio A 

In this section, we discuss FMR by numerically calcu- 
lating the two-dimensional Eq. (|3]) under the condition 
of 87 Rb, namely Co 3> — c 2 > 0. We began calculating the 
spin precessions by applying a n/20 pulse to the ground 
state, whose spins were polarized to the uniform mag- 
netic field H = Hz trapped by V = Muj^x 2 + \ 2 y 2 )/2 
with gfisH/huj x = 20 and an aspect ratio A = u y /Lu x . 

We investigated the dynamics of (F x ) for A = 0.5, 1, 
and 1.5 with and without the MDDI. From t = to 
ir/(20j e H), a 7r/20 pulse was applied. Then, the spins 
were tilted by 7r/20 radians from the z axis with pre- 
cession. After turning off the pulse, the spins precessed 
around the z axis, conserving (F z ). We define the nota- 
tion (Fi)fL x and {Fi)\ = \ a as indicating the expectation 
values of Ft with and without an MDDI in the trap with 
A = A a . 

First, the typical motions of spins are shown in Fig. [T] 
Investigating the time development of (Fi}fL 5 , {F i )fL 11 
and (Fi)fL 15 , we obtained the differences between their 
precession frequencies, as shown in Fig. Q] (a) and (b). 
The differences appeared at frequencies below the Lar- 
mor frequency, given by H. For < t > 2, no devi- 
ation between the precessions was observed, but devia- 
tions clearly appeared as more time elapsed. In order 
to demonstrate that the A dependence was given not by 
H but by i?dd, we show precessions for the same aspect 
ratios without the MDDI in Fig. Q] (c) and (d). The 
precession frequency did not change without the MDDI 
for different values of A. Therefore, the dipolar frequency 
^dd = 7e-Hdd depends upon the shape of the condensate. 

Next, we examined the effects of the MDDI on the 
precessions in Fig. [3J Comparing (F x )f d with {F x )\, 




1 2 30 31 32 



UJ x t bJxt 

FIG. 1: (Color online) The time development of (F x )f d , (a) 
and (b), and (F x )\, (c) and (d). The red solid, blue dashed, 
and green dotted lines show the results of A = 0.5, 1, and 1.5 
respectively. The gray zone represents the duration of a 7r/20 
pulse. 



we observed that the MDDI caused an effective mag- 
netic field, because the frequency of the precession with 
the MDDI deviated from that without the MDDI in 
Fig. M (a) to (f). Assuming that (F,)^ - (Fi)x=i 
is represented approximately to Acos"f e (H + Hdd)t — 
Acos"f e Ht with an amplitude A, we extracted the dipole 
frequency from the waveform. Since the waveform be- 
came — 2Asinwddi/2 sin(o;L +Wdd/2)t, the beat consisted 
of the large frequency ljl + ^dd/2 and the small fre- 
quency Wdd/2. From Fig. [2] (h), we estimated these 
frequencies to obtain ui^d/^L — 6.5, 9,andll x 10 -3 for 
A = 0.5, l,an<il.5 respectively. 

Figure [3] shows the A dependence of uidd/^L- From the 
results, however, we cannot safely conclude that the A 
dependence of the frequencies is given by changing the 
shape of the condensates, since the dipolar frequencies 
may be given by change of the density with the shape. 
FMR in condensed matters has been discussed in con- 
densed matter of uniform density, even with changing 
shape. On the other hand, atomic BECs have tunable 
density and shape. Therefore, our calculations indicate 
characteristic of FMR in atomic cold gases. 



B. Spin current 

We observed spin currents driven by spin diffusion, 
which was caused by a r dependence of the dipolar field. 
Figure H] shows the projections of F onto the x - y plane 
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FIG. 2: (Color online) Comparing the precession with and 
without the MDDI. (a) and (b), (c) and (d), and (e) and (f) 
show the precession for A = 0.5, 1, and 1.5, respectively. The 
solid and dashed lines are {F x )\ and (F x )x ■ (g) and (h) repre- 
sent ((F x }fL , 5 - (F x ) x=0 . 5 )/H (solid), ((F w )? =1 - (F x ) x=1 )/H 



(dot), and {{F X )H 



(F x )\ = i.s)/h (dashed), respectively. 



for A = 1.5 and w x t — 12.7. The precession with the 
MDDI lost homogeneity of the spin directions, whereas 
the precession without the MDDI maintained this ho- 
mogeneity. This is because the precession frequency 
has an r dependence, specifically, u>(r) = j e H e g(r) — 
le (H + H dd (r)). 

The dipole interaction drives the spin diffusion, which 
is shown in Fig. [SJ The figure shows F x /\F xy \ = cos<j) as 
a function of x at y = 0, where <f> is the angle between 
the spin vector and the x axis. In the dynamics with 



Udd 
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FIG. 3: A dependence of Wdd/^L- 




FIG. 4: (Color online) Projection of F onto the x - y plane 
for A = 1.5 and ui x t — 12.7. The figures show the results (a) 
with MDDI and (b) without that. The vectors are nondimen- 
sionalized. 



the dipole interaction for A = 1.5 (a) and 1 (b), the 
spin densities lost their angular coherence, whereas the 
dynamics without the dipole interactions maintained this 
coherence ((c) and (d)). 

The spin diffusion drives the spin current 3k in Eq. 
(1101) . which is shown in Fig. [B] In order to explain how 
the spin current is driven by the spin diffusion, we con- 
sidered the amplitudes of the wave functions tpj = fje lVi 
as, 

^(r,t) = ^y\ l + cosflM))e^ r ^, 



V>o( r >*) 
^-i(r,i) 



n(r,t) 



sin 6(r, t)e 



iip (r,t) 



^^(l-cos^))^- 1 ^, 



(16a) 
(16b) 

(16c) 
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(a) A = 1.5 



(b)A = i 




FIG. 5: Dynamics of a cross-section of F x /\F xy \ at y = 0, 
where \F xy \ = y/F x + F£. From the relation F x = \F xy \ cos</>, 
the parameter represents cos <f>. The results with the MDDI 
(a) and without it (b) are shown for A = 1.5, and (c) and (d) 
show results for A = 1. The x axis are nondimensionalized by 
yfh/Muj x 



and, 
J, 



(1 + COS0) COs(yi — tp )\7tfi 



nh 2 f 
4M\ Sil 

sin 0(1 — cos0) cos((y9_i — tpo)Vtp-i 

2 sin (cos tp r cos tp — cos sin y? r sin y>) \7tp 



2 (cos <p r sin tp + cos sin tp T cos 93) V0 > , (17a) 



- I sin 0(1 + cos 0) sin(<^i - <p ) V^i 

sin 0(1 — cos0) sin(<p_i — ipo)W<p-x 

2 sin (cos tp t sin 93 + cos sin cos tp) Vtp 

2 (cos y? r cos tp — cos sin tp r sin 93) V0 > , (17b) 



nn . , 
4M^ 



COS0) 2 V(y9i - (1 - cos« 



(17c) 



where tp r — (tp\ + tp~i — 2tpo)/2 and tp = (tp\ — tp-\)/2 are 
relative phases. Since the relation tp r = was satisfied in 
our calculations, we used the relation in Eqs. (1171) . and 
the spin density vector formed an azimuthal angle tp with 
the x axis. Then, we derived the spin components F x — 
nh cos <psm6, F y = nh cos tp sin 0, and F z — nh cos 0. We 
can therefore rewrite the spin density currents, 



where the forms show the ground state of the ferromag- 
netic state 23]. The amplitude is represented by n and 
the angle between the spin and the z axis. We intro- 
duced this representation to demonstrate that the spin 
current is derived from the spin diffusion. Of course, we 
confirmed the validity of the ferromagnetic representa- 
tion under the pulse and magnetic field by calculating 
directly. Therefore, it can be utilized for the polarized 
spin state studied in our work. The amplitudes /±i were 
formed to represent F z = nh cos 0, and /o was deter- 
mined to satisfy the relation n = . I^'l 2 - F° r example, 
(n%, no, n_i) = (n, 0,0) led to F z — nh with = 0, and 
(711, no, n_ 1) = (n/A, n/2, n/A) resulted in F z = with 
= 7r/2. The wave function can only express the fer- 
romagnetic states, i.e. the form cannot represent the 
antiferromagnetic state (ni, no, = (n/2, 0, n/2) or 
the polar state (ni,no,n_i) = (0,n, 0). This restriction 
of the wave function is caused by the first representation 
F z = nh cos 0. 

By introducing this representation into Eqs. ^ and 
(fTTj|) . we can redefine as follows, 



F x — nh sin (cos tp T cos tp — cos sin tp t sin tp) , 
F y = — nhsmO (cos tp T sin tp + cos sin ^ cos;/?) 



— (4 cos tp sin 6V<po 

+2 cos tp sm 9 cos 9W tp — 2sin<^V0), (18a) 
nh 2 



4M 



(4sin93sin0V(^o 



J z = 



+2 sin sin cos 0V 95 + 2cos(/?V0), (18b) 

Thfv 

— {4cos0V(p o +4(1 + cos 2 0)V^}, (18c) 
4iW 



which are driven by the gradients of the angles, tp and 0, 
and the phase tp$. In the precessions with MDDI, the gra- 
dients occurred because of the dipolar fields Hdd(r). As 
a result, the spin currents were clearly driven, as shown 
in Fig. |51 For uj x t — 0.12, the spin vectors were coher- 
ent just after the applied 7r/20 pulse (Fig. [6] (a)). The 
spin densities, F x and F y , then flowed to the center of 
the condensates from Fig. [5] (b) to (c) . Then, the den- 
sities reversed, and diffused outward from Fig. \§\ (d) to 
(e). This oscillation was repeated. Of course, we cannot 
obtain the spin current without the dipolar interactions, 
since the gradients of and tp were not caused; the dy- 
namics are shown in Fig. [7] 

In order to investigate the spin fluid dynamics, we cal- 
culated the spin current 3 X for Eq. (1171) . as shown in 
Figs. [5] and M These figures represent 3 X from the pre- 
vious calculations with A = 1 and 1.5 respectively. De- 
spite the difference in the ratio, we observed two com- 
mon properties in these figures. The direction of the 
currents changed rapidly, corresponding to the large pre- 
cession frequency, and the magnitudes changed slowly 
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FIG. 6: (Color online) Dynamics of F projected onto the x 
y plane for A = 1.5 with dipolar interaction. 



FIG. 7: (Color online) Dynamics of F projected onto the x 
y plane for A = 1.5 without dipolar interaction. 



with the small dipolar frequency, as shown in Fig. ITDl 
which shows the time development of the x component 
of 3 x (x = 4,y = 0). This figure indicates that the oscil- 
lation of the current direction occurred with the preces- 
sion frequency, which varied in magnitude with changing 
dipolar frequency. Eq. (ITU) also indicates that the spin 
density was not conserved because of the effective mag- 
netic field. Therefore, the spin currents can be driven 
from a source and sink in the center of the condensates, 
as in Figs. [5] and |H1 The two common properties were 
insensitive to the value of A. However, the change in spin 
density for A = 1.5 exhibited quadratic pole motion in 
a scissors- like mode for mass density {24], which can be 
understood as an oscillation between the spin density mi- 
grating to the y axis from the x axis and back again, as 
shown in Figs. [5] (a) to (c). Therefore, the spin collec- 
tive mode was caused by spin diffusions induced by the 
MDDI. Therefore, the spin current causes the dynam- 
ics of spin scissors-like mode, which was observed as a 
shrinking and expansion of the spin density in Fig. |6] 
The shrinking and expansion were common features for 
A = 1 and 1.5. However, the spin currents were affected 
by the symmetry of the traps, as shown in Figs. [8] and GO 
From the calculations, we expected that the spin cur- 
rent would be observable when using the spinor BECs. 
Recently, spin current is focused from fields of spintron- 
ics. However, it is difficult to observe the spin current in 
metals and condensed matter. Atomic BECs, a macro- 
scopic quantum phenomenon, can show the spin current 
clearly and directly in the dynamics of the spinor den- 
sities. Therefore, we should attempt to observe various 
spin currents utilizing tunable experimental parameters, 
i.e. interaction parameters, trap frequencies, and the 
number of particles. 



V. CONCLUSION 

We investigated the properties of magnetic resonance 
in spinor dipolar BECs by calculating the GP equations, 
obtaining Kittel-like equations as the equations of motion 
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FIG. 8: (Color online) Dynamics of the spin currents pro- 
jected onto the x — y plane for A = 1 with dipolar interaction. 
The vectors are nondimensionalized. 



for the spin density vector. The equations revealed two 
properties. One is the dynamics of the spin fluid, and 
the other is precession under the effective magnetic field 
consisting of the external magnetic fields and the dipolar 
fields. The magnetic resonance with the properties of the 
spin fluid was characteristic of this system. 

In order to extract properties from the GP equations, 
we studied the law of conservation of spin density cur- 
rent without effective magnetic fields by first deriving 
the continuity equations from the GP equation, obtain- 
ing representations of the spin current. Second, we ana- 
lytically evaluated the precession dynamics described by 
the Kittel equations derived from the GP equations us- 
ing a single-mode approximation, where the Kittel equa- 
tions show conventional FMR. The analysis clearly in- 
dicated that the origin of the FMR in the BECs is like 
the dipolar field, whereas the origin of the resonance in 
the Kittel equations for condensed matter is the demag- 
netizing field. Comparing the FMR of the BEC with 
that of the condensed matter, we concluded that the ori- 
gin of the resonance was not the spin exchange interac- 
tion that causes magnetism in condensed matter, but the 
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FIG. 9: (Color online) Dynamics of the spin currents pro- 
jected onto the x—y plane for A = 1.5 with dipolar interaction. 



anisotropy of the MDDI. Finally, we numerically calcu- 
lated the GP equations, representing the dynamics with 
the two common properties. The characteristic dynamics 
showed that the effective magnetic field introduced spin 
diffusion into the Larmor precession, driving the spin- 
current-like scissors modes. 

The relation between the spin current and FMR has 
not yet been discussed for typical FMR. Therefore, it is 
important to study spin current in condensates. We also 
believe that the study of spin current will be useful for 
the development of spintronics, because it is difficult to 
directly observe spin currents in condensed matter spin- 
tronics. 
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FIG. 10: Dynamics of the x component of at x — 4 and 
y = 0. The inter figures are the results for ui x t = to 4. 
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